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ABSTRACT 



Aims. Dissipative phenomena occurring during the orbital evolution of a dwarf satellite galaxy around a host galaxy 
may leave signatures in the star formation activity and signatures in the colour magnitude diagram of the galaxy stellar 
content. Our goal is to reach a simple and qualitative description of the these complicated phenomena. 
Methods. We develop an analytical and numerical technique able to study ram pressure, Kelvin-Helmholtz instability, 
Ray leigh- Taylor and tidal forces acting on the star formation processes in molecular clouds. We consider it together 
with synthetic colour magnitude diagrams techniques. 

Results. We developed a method to investigate the connections existing between gas consumption processes and star 
formation processes in the context of the two extended-body interaction with special attention to the dwarf galaxies 
dynamical regime. 



Key words, tidal forces, ram pressure, Ray leigh- Taylor, Kelvin-Helmholtz, 
mation processes, stellar populations, colour magnitude diagrams 
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1. Introduction 

We investigate the role of dissipative phenomena in the star 
formation history of a dwarf satellite galaxy orbiting around 
a primary galaxy. Our purpose is to quantify in a simple 
way the connections between the gas consumption processes 
(e.g., ram pressure, Ray leigh- Taylor and Kelvin-Helmholtz, 
tidal forces) and the star formation processes in the context 
of the interaction between two extended-bodies with special 
attention to the dynamical regime of dwarf galaxies. 

Dwarf galaxies offer a beautiful example of interac- 
tions between a m ajor galaxy and its satellites (e.g., 
iPccblcs 199 J, |1980D . While orbiting these satellites are 
affected by environmental effects that shape their evo- 
lution and may a ffect their star formation history (e.g., 
iGrebel et al.ll200"3[) . Their star formation history can be in- 
ferred from their colour magnitude diagrams (CMD), see, 
e.g. the recent analysis by .Orb an et al^ (2008.) of Hubble 
Space Telescope data for the galaxies of the Local Group 
(LG). Generally, environmental effects are important in 
groups and clusters of galaxies and are believed to con- 
tribut e to the origin of t he density-mor phology r e lation 
(e.g., lOemled Il974t iDressler 1980: Dressier et al.l 



IValluril Il993t iBvrd fc ValtonenI |1990D thus induc i ng star 
form ation episodes (e.s.. iHenriksen fc BvrdI 119961 : lEvrardl 
Il991h . Tidal compression is also responsible for the so-called 
gala xy-harassment phe n omen o n in t he case of dwarf galax- 
ies ("Moore et al.l |l996l Il998l I1999D or late- type systems 
(e.g., Gncdin 20of~ In t his context, three- dimensional nu- 
merical simulations fe.g.. lQuilis et al.]|2000() show that ram 
pressure can be effective in removing up to 100% of the 
atomic hydrogen from luminous galaxies within 10^ yr. 

With respect to the LG dwarf galaxies, the role of 
the external mechanisms of gas removal has been inves- 
tigated by means of different ap proaches and with con- 
flicting results (for a review see iGrebel et al] 12003") and 
in relation with the MW we cannot neglect the dissi- 
pative role of tidal forces by the deep MW potential 
and its hot coronal plasma. A large amount (~ 5 x 
lO^^AffT)) of hot gas is thought to exist in an extended 



(~ 200kpc) hot diffuse halo (see 


e.g., 


Briins et al.l 120001: 


Sembach et al.ll2003l: Putman et al.ll2003 


ilMaller fc Bulfockl 


20041: Fukueita fc PeeblesI 20061 


Kaui 


"mann et al.l 20061: 


Sommer-LarsenI 2006f). Ram oressure can both enhance 



1997) 



or to the Butcher- Oemler effect (iButcher fc Oemleij 11978 ; 
ICouch et all Il998l : iGoto et al.l |2003D. The mechanism of 
tidal compression exerted by the gravitational potential of 
a galaxy cluster is known to increase the velocity disper- 
sion of the molecular clouds falling in disk galaxies (e.g.. 
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star for mation by compression of the interstellar medium 
(ISM) (jPuiita fc Nagashimalll999D and reduce it by strip- 
ping the gas from the galaxy ( Gunn fc Gottl[l972n . Also, 
tidal forces can trigger episodes of star formation. These 
are the major external competitive effects that we want 
to investigate in the context of 'dwarf galaxy-massive 
galaxy' interaction. In addition to this, there are also sev- 
eral internal mechanisms that govern the gas consump- 
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tion, such as gas expulsion by star formation processes 
Burkert & Ruiz-Lapuentc '1997| iFerrara &: Tolstovl 
Mac Low fc Fcrrara 1998) giant HII regions (e.g. 



2000: 



Davies fc PhilliDPsill988l: lElmegreen fc Hunted f2000l ) . The 
role of the hot gaseous halo surrounding the MW has of- 
ten been studied in relation with the Magellanic Stream, 
a long HI filament leading and trailing the Magellanic 
Clouds (see e.g., .Murah 2000; Nidever et al. 2010. ), where 
the authors constrained the density of the coronal gas to 
be n < 10~^cm~^ at the distance of 50 kpc in order to 
let the Stream clouds survive up to 2.5 Gyr or at least a 
minimum of 0.5 Gyr. IStanimirovic et all ()2002f ) suggested 
a limit of n = 3 • 10~'*cm~^ at the same distance, which 
is essenti ally in agreement wi t h the more recent measure- 
ments by 'Heitsch & Putman' ('20091) . Previous determina- 
tions by Kalbcrla & Kcrp (1998) based on the ROSAT X- 
ray observations agreed on the same limit for the MW halo, 
whereas mode ls predicting dwarf galaxies to lose m ost of 
their gas fe.g.. lGallart et al.ll200ll:lEinasto et~aIlll974D need 
a gas density one order of magnit ude higher. On the other 
hand, iGrcevich fc PutmanI ()2009l) argue that ram pressure 
is one of the primary causes for the gas deficiency of dwarf 
sate llites within 270 kpc around the MW a nd M31 (see 
also Ivan den Bergh|[T99l iGrebel et"all I2003D . Finally, we 
assume that the present-day density of the intergalactic 
medium beyond the LG halos is too low for ram pressure 
to be effective fe.g.. iQuilis fc Moore|[200l . 

The subject of our analysis will be an arbitrary ideal- 
ized dwarf galaxy orbiting around the Milky Way (MW) 
galaxy. The differential equations involved in our study are 
generally coupled in a non-linear way thus not allowing us 
to consider the global evolution of the processes as a linear 
superposition of the single effects considered alone and then 
summed up. Nevertheless, for the sake of clarity, it is cus- 
tomary to subdivide the effects involved in the evolution of 
a dwarf galaxy in two categories: 'internal' and 'external'; 
and then to consider their superposition. We will follow this 
tradition also in this paper. Among the internal effects, stel- 
lar evolution processes are the most important ones because 
they ultimately govern the location of the stars in a CMD. 
The evolution of the stars and their feedback is mainly re- 
lated to their mass and chemical composition. In our ap- 
proach, we will carefully consider stellar feedback on the 
dwarf galaxy stellar populations in dependence of the mass 
of the stars (see Section [2. ip . Finally, in our approximation 
we want to consider several processes of external nature 
as well, mainly tidal compression, ram pressure, Rayleigh- 
Taylor and Kelvin-Helmholtz instabilities. The reasons for 
focusing on these effects are the following. We are mostly 
interested in the effects of the orbits on the generation of 
observable CMD. To approach this problem we will work 
out a spherical approximation for the shape of the ideal- 
ized dwarf galaxy considered as a whole. On the one hand, 
this strong approximation will make us lose track of the 
orbital evolution of the stars inside the galaxy itself, i.e. 
we will not solve the equation of motion (EoM) for in- 
dividual stars inside the dwarf galaxy (see Section 12. 3p . 
but this approach will allow us to ne glect the effects of 
tidal stripping and tidal stirring (e.g . jMayer et al. 200lt 



i g ana tiaai stirrm g (e.g ., 'Mayer et ai.llzuuit 
1 120031 iRead et all 120061: Klimcnt ~ki et al.l 



Pasetto et al.l 

2009t) giving more agility to our code (sec the code per- 
formance in Appendix A ) as compare d with standard N- 
body techniques (see, e.g.. lMaveill2010l for a recent review). 
Finally, in order to simplify further the orbit integration. 



we will limit ourself to non-penetrating encounters with the 
MW, i.e., we will limit our analysis to the dwarf galaxy orbit 
pericenter, Vp, greater than ||rj,|| > 20 kpc, thus neglecting 
partially the problem of the indetermination of the parame- 
ters of the MW potential, the role of the dynamical friction 
and role of tidal shocks due to the passages trough the MW 
disk (which are not the target of the present paper). Being 
mainly focused on the LG dwarf galaxy orbits, we will not 
consider galaxy harassment because it is muc h less efficient 
in low density groups such as our LG (e.g., iMoore et al.l 
[T996h . 

The paper is organized as follow: In Section|2]we present 
our model for gas the consumption and star formation in 
dwarf galaxies taking into account both internal and en- 
vironmental effects that would affect the star formation 
history while the dwarf galaxy is orbiting about a host- 
ing galaxy (the MW in this paper). In Section [3] we shortly 
summarize the key lines of the population synthesis tech- 
nique to generate colour-magnitude diagrams out of which 
the history of star formation is inferred. In Section |4] the 
technique is applied to a model of a dwarf galaxy in which 
the star formation efficiency is altered by the environmental 
effects during the orbital motion and from which a synthetic 
CMD is generated and shortly analysed. A brief discussion 
of the results and some general conclusions are presented 
in Section [5] 

Finally, we made extensive use of the Appendix to 
present technical test on the code or to explain the ana- 
lytical aspects of the work that can be postponed at a first 
reading of the paper. 



2. A model for gas consumption and star formation 
in dwarf galaxies 

2.1. Molecular clouds evolution and non gravitational heating 

To link the orbital motion of a satellite about its host 
galaxy to the environmental effects and the star forma- 
tion efficiency in the satellite itself, we must set up a fast 
algorithm allowing us to explore a large space of initial 
conditions for the orbital motions. This immediately rules 
out a multi-phase fully hydrodynamic treatmen t of the gas 
fe.g.. lOuilis et al.l l200dl; iTonnesen et al.l l2007t) because it 
would be too time consuming in particular if a large vol- 
ume of initial conditions (i.e.) of the orbital motions has 
to be explored. Therefore, we prefer to follow the m ethod 
proposed bv lFuiital ^9M) and lPuiita fc Nagashimal (|1999D 
from whom we take the model of star formation to be in- 
tegrated by extensive calculations of orbital motions and 
synthetic CMDs of the stellar content of the dwarf galaxy. 

In brief, stars are known to form in gas clouds (ei- 
ther atomic or molecular) whose mass distribution is un- 
certai n but believed to span at l east a factor ^ 10^ in 
mass (jElmegreen fc EfremovllT997l ). We group the molec- 
ular clouds according to their initial mass Mi\Mi g 
[Mniin, Minax] and assume M^in = IO^Mq and M^ax = 
1O^M0. Since this mass interval is very large, we split it 
logarithmically and resolve the mass intervals with a ra- 
tio a between any two subsequent intervals, i.e. we solve 
the recurrence relation logj^g ^^jf^ = a, with the boundary 

condition AIq = Mmin for i=0, as Mi = 10"*Mmin (a is a 
rational number). The numerical convergence is confirmed 
for a = (see Appendix A for details). 
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We then define for simplicity a "mass class" within the 
generic interval Mi = Af^+i — Mi or its fraction over a total 

mass Aftot as fi = ^^Mtot' assume that the evolution 

of each mass class Mi representative of each mass interval 
is governed by the fraction of gas that is ejected from the 
stars, say their gas ejection rate i?star and by the amount 
of recycled molecular gas with recycling rate i?moi- 

To determine i?star we remember that short lived-lived 
and long-lived stars contribute to i?star in a different fash- 
ion. If r (m) is the amount of mass returned by a star of 

mass m, fi (m) — ^-^^ is the fractionary return mass. Then 

[SI 

the gas ejection rate due to stars with short lifetime, iigtar? 
will depend on the mass distribution of the stars, i.e. the 
time independent initial mass function l — l. (m) expressed 
as a mass fraction, and the star formation rate ip at the 
moment in which the star is born 'tp (t — i (m)^ with i (m) 
being the lifetime of a star of mass m. Supposing that the 
typical age of a galaxy is, for instance, to — 12 Gyr and 
looking for mass of stars whose lifetime is small compared 

rg] 

to the age of the galaxy (say 1/10 to), RitL given by 



[S] 



(t) = / ip (t — i {m)'j fj, (m) L (m) dm, (1) 



where miow = 2.3M0 is the star with t{m) ~ O.lic and 
niup is the most massive star we want to consider 



IOOMq. All the less massive stars contribute to the long- 
lived part of i?star indicated as R^^l^,-- The contribution of 



these stars can be written as R, 



[L] 



- (0). ^0 is 
an artificially fixed star formation rate for the galaxy we 
are considering and we consider the lifetime corresponding 
to the transition mass from long-lived to short lived stars 
ttr = O.Ug Gyr. Recollecting our previous terms, we get 



-fl-star — ^star 



R. 



[L] 



Cr [V' (i - i) - V' i-t)] fiim)Lim)dm, 



(2) 



where 9 (x) is a step function {6 — for x < 0, 6* = 1 

otherwise), and i?star is indeed not time dependent for t > 
0. 

Furthermore, we define the molecular recycling rate, 
i?inoi, as the complement of the star formation efficiency, 
1 — e where e is a general function of the molecular cloud 
mass M and the pressure P, i.e. e = e{M,P). The star 
formation efficiency in turn can be reduced by the presence 
of an external gas consumption mechanism, e.g., HI phase, 
ram pressure, Rayleigh- Taylor instability, tidal compression 
etc. (see the Section [2.41 below for a detailed discussion). 
Therefore, for each mass class, i, we can define 



i?mol (i) =^(l-e(Af„F)); 



M, (t) 

■m,p) 



(3) 



which depends on all the molecular clouds i, and where 
T {Mi, P) in the destruction time of a molecular cloud with 
mass Mi and pressure P (see below for an extended dis- 
cussion, and Appendix B for the analytical formulation). 
In principle, it is now simple to introduce in Eqn. a re- 
ducing factor in i?moi in order to lower the fraction 1 — ScS 



when the stronger energy injection by supernovae, photo- 
ionization of molecular clouds, or the effects of strong stel- 
lar winds need to be investigated. We have chosen to focus 
the present work mostly on external effects on the gas con- 
sumption or gas removal. 

Finally, the time variation of Adi simply depends on the 
life-time r of the molecular cloud according to the system 
of integro-differential equation^B: 



dMi 
dt 



fi (-Rstar + -Rmol) 



T 



(4) 



which relates the rate of change of each mass class to 
all the other molecular classes though the gas ejection rate 
of the stars i?star, and the recycle rate of molecular clouds 
R-mo\ (where we again use the step function defined after 
Eqn. ^ to take into account that the mass is a positive 
definite quantity). 

Once the system of equations ^ is considered together 
with Eqns. (HJ, ^ and we can estimate the star for- 
mation rate (SFR) as 



V'(t)=5^£(M„P) 



M, {t) 

■{M^,py 



(5) 



■th 



all the other — 1 molecular clouds as seen in Eqn ([TJ and 
represents the target of our computation. 

In the previous system of equations the only functions 
that needs to be adopted from the literature are the star 
formation efficiency of a molecular cloud e {Mi, P) and the 
destruction time r {Mi, P) of a molecular cloud with mass 
Mi and pressure P. Then, once these function are assumed, 
the system can be integrated once ipo is assumed on the ba- 
sis of the system we want to analyse (e.g., globular clusters, 
dwarf galaxies, spiral galaxies etc.). 

Star formation most likely responds to the pressure vari- 
ation in the local interstellar medium, and in our case 
it is worth recalling that the pressure P = P{x,t) de- 
pends on time and relative position with respect to the 
host galaxy (with the barycentre centred in the origin O 
of the inertial reference system 5*0 with vector radius x) 
and the environment. The relation between star formation 
efficiency, e, mass M and pressure P of molecular c louds 
has been studied e.g . , inlEl megrec n &: Efremovl (|1997t ) and 
iKrumholz fc McKed (|2005l ) showing how a pressure in- 
crease causes an increment of e and a subsequent decrease 
of the lifetime of a cloud r. These relations are independent 
of the physical origin of the pressure and can be safely ap- 
plied to our treatment of the SFH during the orbital motion 
of a dwarf galaxy. We simply assume these relations from 
the literature. A plot of their values can be seen in Fig. 
[1] for the range of values adopted in the case of the dwarf 
orbiting around the MW. We refer the interested reader to 
further details and their analytical definition in Appendix 
B. 

Finally, we mention that the dynamical response to 
gas removal by SN explosions in dwarf galaxies has been 
repeatedly invoked as a mechanism to suppre ss the star for- 
mation in these objects (e.g., So merville fc Primack 19991 : 



^ To reduce the degrees of complexity of the equations, we will 
not consider cosmological effects or the explicit dependence of 
the equations on the redshift. 
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0.0 0.1 0.2 0.3 0.4 

t{M,P) 

Fig. 1. Star formation efficiency vs. fifetime relation. The 
surface is an interpolation with 20 meshes of constant pres- 
sure run horizontally and linearly for P € [0.001, 1] Pq from 
top to bottom. Ten mesh lines runs vertically for constant 
mass, logarithmically for Mi € [lO^, 10^] Mq from left to 
right. 



Cole et al.ll2000HDe Lucia et al.l[200l iNaeashima &: Yoshiil 
2004 iKang et all 1200511 These results are never- 
theless in contrast w ith observational p ro pertie s 
as pointed out in Steiiimetz fc Navarro! (Il999l). 



iNavarro fc Steinmet j (l2000af). Navarro fc Steinmetzj 



(|2000bD . iPortinari fc Sommer-LarsenI (|2007t l suggesting 
that numerical recipes can be improved as, e.g., recently 
shown by Pasetto et al. (2010) using prescriptions based 
on chemo-dynamical studies. These authors made evident 
that the gas survival in spherical dwarf galaxies is perfectly 
compatible with the expected SN rate over the whole 
mass range of interest for dwarf galaxies (see their Table 
1). This reasonably stems from the reduced SN feedback 
expected in the local environment {Cioffi & Shull IQOlh 
differ ent from dwarf to spiral galaxies (Bradamante et aD 
|1998( ) and it explains in a natural way the co-existence of 
old stellar popu lations a nd a reservoir of neutral gas (e.g., 
iPhiUipDS et"al]ir99Q:iPasetto et al.ll2010t) . 



2.2. Galaxy model and plasma distribution 



The most 
ing of the 
physics 



important contributions to the understand- 
plasma corona of the MW come from par- 
Muralil l2000t ) and cosmology (e.; 



tide physics (e.g., . . v- o , 

iFukueita fc Peebles! l2006i r Nevertheless, the literature is 
rich in attempts to model the distributions of the hot 
gas density (which because of the high degree of ion- 
ization refers uniquely to electron density) in differ- 
ent environments e.g. in the LG (R asmussen & Pedersen 



2001 



19991) 



iBinnev fc EvansI l200ll iMalonev fc Bland-HawthornI 
by considering different constraints, in clusters 
of galaxies by fitting the X-ray observation s with /3 - 
models (e.g.. ICavaliere fc Fusco-Femianol 19761) and by fit- 
ting the MW hot coronal gas (e.g., Moore fc David 1994 
Fukugita & Peebles 2006; Ric hter et al.ll2008HSavage et al.l 
2003.: Tumlinson fc Fang ,20051 ). In a hierarchical galaxy 
formation scenario, after the initial collapse, the con- 
tinuing accretion of gas-rich fragments can produce a 
diffuse hot gas halo that surroun ds the galaxie s and 
fills the dark matter potential (| White fc ReesI Il978t 



I White fc Frenk 1991) not necessari ly with the same radial 
profiles (jFukugita fc Peebles! |2006[ ). This scenario agrees 
also with the case of clusters of galaxies where the X- 
ray distribution of the intra-cluster medium is signifi- 
cantly flatter than the dark matter dist ribution, in partic- 
ular for low-temperature clusters (e.g., lArnaud fc EvrardI 



5 X 10" 



2x10" 



50 

R [kpc] 



Fig. 2. Electron number density for our coronal model. For 
comparison the dashed line show s the electron number den- 
sity model taken from literature (jFukugita fc Peeblesll2006l 
their 'model 2'). 



119991: lEttori fc Fabian|[T999t iMohr et al.lll999[ ). due to non- 
gravitational heating. 

In our modelling, we assume that the density distribu- 
tion of this gas settles down in hydrostatic equilibrium un- 
der the effect of the gravity of the dark matter component 
of the host galaxy (the MW in our case) with a temperature 
around lO^^'^K (or ~ 0.086-^0.86 keV as common in plasma 
physics). The ion component provides the density and it is 
sel f-consistently added to the galactic potential described 
bv lPasetto et al.l (j2011) to which the re ader should refer for 
a the d etailed description. Followin g iFukugitafc Peeblei 
(|2006D (jAnderson fc BregmanI [20Tol see also,), we adopt 
Y = 0.25 for the primordial abundance of helium, with a 
sound speed Vs — y^jZkBTe/rrii = 600 km s~^ (with Tg as 
the temperature in eV, rrij as the ion mass, Z as the charge 
state, ks as the Boltzmann constant and 7 as the adiabatic 
index (see also iKippenhahn fc WeigertI 119941) . The result- 
ing model for the electron density is shown in Fig [51 Other 
approaches based on polytropic equations of state cannot 
add any significant insight to the problem due to the lack 
of observational constraints on the plasma temperature's 
spatial gradient or rotation in the range of distances we are 
interested to sample. 



2.3. Tidal approximation for a pressure equation in nan 
inertial reference systems. 

We limit ourselves to the simple case of a spherical dwarf 
galaxy whose size is given by a typical scale radius, r^, e.g., 
the tidal radius or the effective radius, or any other ra- 
dial scale length that might be obtained from observational 
data. If the dwarf galaxy is not a sphere, then the direction 
of motion must also be specified in its orbit around a mas- 
sive companion in relation to the shape of the dwarf galaxy, 
and more sophisticated formulations are required. 

We want to develop a formalism for the pressure acting 
on a dwarf galaxy while orbiting through the hot coronal 
gas of the MW (Section l2.2p . Quantifying the pressure is 
necessary in order to investigate the effects on the star for- 
mation rate as explained in Section [2.11 

The pressure can be obtained from the Navier- Stokes 
equations in a potential-flow approximation once these 
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equation are written in a non-inertial reference system, 5*1 , 
comoving with the dwarf galaxy. This approach is based on 
the classical assumption that the interstellar medium (ISM) 
of a dwarf galaxy is much denser than the hot intergalac- 
tic medium (HIGM) of the MW. If we consider that the 
HIGM has a temperature greater than 10/100 times the 
dwarf galaxy ISM temperature (see Section 12. 2p then we 
can neglect the effect of the dwarf potential on the HIGM. 
Moreover the ISM will be considerably denser than the 
HIGM so that, to a first approximation, th e flow of the 
hot gas is like that around a rigid body fe.g.. lNulsenlll982l 
and references therein). In this case, we can proceed by 
adopting the standard formulation for the potential flow 



(6) 



(e.g.. [Chorin k Marsden|[l979l:lBatcheloijl2000t) . Here (a, b) 
is the inner product between two arbitrary vectors a and 
b, the subscript Vq refers to the fluid velocity of the po- 
tential flow given by Vq = V(p„Q, v is the dwarf galaxy 
velocity in an inertial reference system 5*0 centred on the 
MW, and ^ the vector radius in the non-inertial reference 
system Si whose origin is in the barycentre of the dwarf 
galaxy O' . Hereafter the prime-notation refers to the non- 
inertial reference system Si. The integ ration of the Euler 
equation is a standar d procedure fe.g.. iLandau fc Lifshitzl 
I1959I : lBatchelor1l200G( ) that leads to the classical Bernoulli 
equation that we would need to solve for the pressure P, 
but instead, we write it in a non-inertial reference system 
by adding the extra term {A, ^) that takes into account the 
non-inertial character of the reference system Si (where A 
is the acceleration of a fluid element due to the motion of 
the galaxy). We get: 



dt 



P 
P 



= f{t)-^g-{A,S,) 



(7) 



In this equation, p is the coronal gas density distribution 
and $g the gravitational potential acting on the fluid ele- 
ments. We use a prime to indicate that the potential flow 
refers now to a non-inertial reference frame Si as well as 
its velocity Vq. This equation closely resembles the clas- 
sical Bernoulli equation but fo r the presence of the n on- 
inertial term (see Section 3.1 of iPasetto fc Chio"sill2009l for 
a more extended discussion). Finally, we fix the arbitrary 
/ {t) (constant over all the streamlines for inviscid irrota- 
tional fluids) by imposing the boundary condition that the 
hot gas corona is in hydrostatic equilibrium. We easily ob- 
tain 



fit) 



(8) 



The last approximation we want to introduce before 
solving Eqn. ([7]) for the pressure P is the tidal approxi- 
mation. Because we are interested in the pressure acting in 
the regions close to the dwarf galaxy, we perform a Taylor 
expansion of the total gravitational potential $g (a;) acting 
on the generic element of the fluid located at x close to the 
barycentre (located at ilp) of the dwarf galaxy: 



$g {X) ~ $g (i?0) + V$g (Ho) (3; - i?0) + 



(9) 



so that X Rq — 0$, with O G SO (3) is the generic rota- 
tion matrix between the orthonormal basis of the inertial 



reference system Sq and 5*1 with transpose , and the 
approximation holds for ^ <^ Rq and ^ x. Considering 
Eqn. © and that A — — V<f>g, we easily obtain from Eqn. 
^ the required equation for the pressure in the tidal ap- 
proximation 



P = p 



y (4 - 9 sin^ d) 



$g iRo)+r^,{0'^TOfi,n) 



(10) 

where = is the square of the norm of the veloc- 
ity vector of the dwarf galaxy, i} is the angle between the 
direction of the velocity vector and the normal h (e.g., 
iLandau fc Lifshit3ll959[ ). and the equation is evaluated at 
the scale radius: ^ = r^n with h = |||^. Finally, we have 

introduced the tidal tensor T = —H{^g (Rq)) defined as 
the negative to the Hessian matrix H in the Taylor expan- 
sion of Eqn. (|9]) stopped at the first order in the potential 
and evaluated at the dwarf galaxy barycentre position. 

Equation ([T]) is completely determined once the orien- 
tation of the reference system with respect to is de- 
fined, i.e. when we define the rotation m atrix O. We use the 
Cayley-Klein parameters (e.g., Goldstei n et al.ll2002f) to ori- 
ent Si with the Frenet-Serret reference frame {eT,eN,^B} 
once the orbit is numerically integrated where the tangent 
vector Et, the normal e^v and the binormal eb are its in- 
stantaneous associated unitary vectors. From the orbit we 
can infer the kinetic pressure P = P {x,t), where x = x (t) 
is the solution of the EoM in the inertial reference system 

50 centred on the host galaxy de fined by the MW gr av- 
itational potential $3 as in, e.g., IPasetto et al.l ()2011[ ) at 
tih = 0, i.e. the present-day potential (with iib look-back 
time) . For simplicity in the following we are neglecting the 
time dependence of the gravitational potential. 

We highlight a few characteristics of Eqn. ([TU)) : 

— It contains an explicit dependence on the size of the sys- 
tem Ts. This parameter can be easily changed from the 
tidal radius to the half-mass radius, or can be related to 
the light profile depending on the specific applications. 

— It shows a dependence on the direction n along which 
we want to estimate the pressure. 

— It reduces to the standard results in the special case of 
the absence of a gravitational field, in the case of an 
inertial system of reference and for circular orbits. 

In Figure [3] we show the pressure in the reference system 

51 for a dwarf galaxy orbiting with an eccentricity e = 0.5 
and pericenter Vp = 200 kpc in the MW external potential. 
This model of the galaxy is later referred to as model E050 
(see further details in Fig. [5]) . Along each direction perpen- 
dicular to an ideal surface enclosing the ISM of the dwarf 
galaxy (indicated by the blue color) a cylinder or piston 
is plotted. Its height is proportional to the intensity of the 
pressure normalized to the value at the stagnation point 
(the ideal point where the flow velocity vanishes). Due to 
the complexity of the terms in equation ()10p the behaviour 
of the cylinder is clearly not linear. 

Finally, we can further refine our model in two way: 

— By considering the variability of the scale radius that 
can occur during the evolution of the dwarf galaxy in its 
tidal interacting with the host system. In this case, due 
to the linearity of the Laplace equation for the potential 
flow, we can simply add the potential flow of a reducing 
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Fig. 3. Plot of the pressure acting on the "surface" of a 
dwarf galaxy (blue sphere) with orbit E050 at the look- 
back time tib = 9 Gyr (see the text for further details and 
Fig. [5]). Orange pistons represent the pressure intensity on 
the dwarf, the higher the piston the higher is the pressure. 



size dwarf 



(11) 



to the previous potential flow of Eqn. ^ . In this case to 
deduce the pressure term is more tedious but straight- 
forward and we omit it. 
— By considering that the orbits of dwarf galaxies are 
different in nature according to whether their velocity 
is subsonic or transonic/supersonic. When the galaxy 
velocity becomes comparable with or exceeds that of 
the sound in the inter galactic medium, effects due to 
the compressibility of the fluid become important and 
shock waves can occur wi thin the g alaxy or backward 
throughout the ISM (e.g.. iLandaufc Lifshitz 1959). In 
the literature these concepts find application in prob- 
lems related to the free fall of galaxies in rich clusters 
of galaxies (e.g.. lHenriksen fc B vrd 199(| their Eqn. 6). 
Nevertheless we will not encounter this situation in the 
dwarf galaxies orbiting through the hot MW coronal 
gas that we are going to analyse since the velocity and 
density involved are far from the supersonic regime. We 
refer the interested reader to Appendix IA.3I for the test 
of the pressure equation and its relation with the stan- 
dard supersonic formalism. 

In the panels of Fig. |4l [5] and [6] we plot three represen- 
tative orbits and the relative pressure computed with Eqn. 
(fTU|) along the direction of motion er- The pressure is plot- 
ted for the unitary scale parameter rs — I, normalized pres- 
sure and time in Gyr. As arbitrary normalization value for 
the pressure we have chosen the standard literature value 



for pressure from the random motion of molecular clouds 
in the solar neighbourhood, Pq = 3 x W^kn cm~^K with 
ks the Boltzman n constant (e.g., lElmegreen fc Efremovl 
119971 : lFuiitalll998l and Appendix B). As expected the max- 
imum pressure experienced by the dwarf is at the point 
of closest passage in the orbit mostly as a result of the 
higher velocity of the dwarf. The orbit selection is based on 
the eccentricity. We assume a starti ng distance of 2 00 kp c 
in the MW potential described by iPasetto et al] (|2011[ ). 
then we search for initial conditions of the orbit in the 
MW plane supposed to be at rest in starting from 
Wq = {xq, Vo} |e — {0.25, 0.50, 0.75} with initial conditions 



xo = {0,200,0} kpc, e 



rp/ra+ rp where r^, rp are 



apo-centre and peri-centre of the orbits respectively. We 
will refer to the 3 orbits evolved in the MW model intro- 
duced in the previous Section by their eccentricity E025, 
E050, E075. 

2.4. External gas consumption effect: ram pressure stripping, 
Kelvin-Helmholtz and Rayleigh-Taylor instability 

Given the prescription for the pressure along an assigned 
orbit, we can reconstruct the SFH of the galaxy by resolving 
the system of equations However, we need also to con- 
sider the effects of external instabilities that can affect the 
gas consumption processes and, as a consequence, the SFH 
of a dwarf galaxy as well. We consider two distinct types of 
stripping processes: (i) instantaneous stripping due to ram 
pressure and Rayleigh-Taylor instabilities and (ii) continu- 
ous stripping due to the Kelvin-Helmholtz instability. 

The approach adopted to deal with the kinds of strip- 
ping is different. For the ram pressure and Rayleigh-Taylor 
instabilities we extend the classical recipe by making use 
of the formalism developed in the previous section. It has 
the advantage of taking the galaxy size and the direction 
relative to the motion into account analytically. In this 
way, we gener alize the literature re sults, see e.g., the clas- 
sical study by Gunn fc Gott' (1972^) and the recent works 
by Rocdigcr & Hciislcr (2005) for disk galaxies or those 
bv iMori fc BurkertI (|2000j) on the thermal pressure in the 
galactic centre. 

Stripping of gas by ram pressure occurs when the pres- 
sure times the unit area, i.e. —Pg^tsfi'dA, is larger than the 
restoring gravitatio nal force for a un i t mas s of gas (HI or 
molecular, see e.g., iKrumholz et al.l (|2009( )) of the dwarf 
galaxy clouds (with $dg, M^g and pdc the potential, mass 
and density respectively in 5*1 of the dwarf galaxy) : 



sir) 



Mdg (e,i) =4^ 



e 



(12) 



For simplicity, the time dependence is only in the total mass 
derived from Eqns. ^ whereas the scale of the structural 
parameters remains unchanged. 

This formalism can be easily applied to study the 
Rayleigh-Taylor instability that in order to occur requires 



e 



Mg.s {i,t) < ^CdP, 



(13) 



where the drag coefScient, Crf, is taken to b e equal to one for 
simphcity fe.g.. iRoediger fc Hensleijr2008D . It is also inter- 
esting to note that the formulation adopted in Eqn. ([TU)) is 
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O.2O1 




71b [Gyr] V-I 

Fig. 4. For a test satellite orbit with eccentricity e=25 and baryonic mass Alt, = 10'^ Mq we show the path (upper left 
panel), the resulting pressure profile (upper right panel) normalized to Pq defined in the text, the star formation rate 
(lower left panel) and the resulting CMD (lower right panel). The colour bar is plotted only in the upper left panel, red 
for the oldest stars at the beginning of the orbital path t = 0, blue for the youngest stars born after 10 Gyr of orbital 
evolution. In the pressure panel, tn, is the look-back time in Gyr. In the CMD the size of the diamonds is proportio nal to 
the m etallicity Z of the stars as obtained the Padova CMD simulator YZVAR fed with the stellar models of .Bertelli et al.l 
(|2009t l (metal rich stars are indicated by large diamonds). 



"local" , in the sense that it contains information regarding 
the direction n of the point where the condition Eqn. 
holds. In this way we can take into account of the direction 
of observation fi. We see here in analytical way, which dwarf 
galaxy parameters are relevant for the Ray leigh- Taylor in- 
stability. This effect has a dependence on the velocity of the 
dwarf the scale length and the direction of the motion 
trough the function sin z9. As soon as we look away from the 
direction of the velocity vector (i.e. ■(? increases) there are 
zones of the dwarf galaxy that are more stable against this 
type of instability i.e., the pressure in Eqn. (jl3p decreases. 



At the same way we see that the more massive the galaxy 
is, the less efficient is Ray leigh- Taylor stripping. 

Finally, the gas surviving the instantaneous ram pres- 
sure stripping can be later removed by the Kelvin- 
Helmholtz instabilities occurring at the interfaces between 
the ISM in the dwarf galaxy and the HIGM, even if their 
effects are partially suppressed by gravity. Since we are not 
inte rested in wakes we can saf ely neglect the ICM viscos- 
ity (jRoediger fc Briiggenll2008h and consider this stripping 
process as continuous and mostly affect ing the gas in the 
flow past the dwarf galaxy ()Nulsenlll982f) . This effect is ex- 
pected to decrease the total mass of the dwarf galaxy in the 
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Fig. 5. As in Fig. Ubut for an orbital eccentricity e=0.5 of the satellites. The pressure panel refers to the stagnation 
point, i/f, is the look-back time in Gyr. For tuj = 9 Gyr the surface pressure is represented in Fig. |3l 



Eqns. (U) according to 



dM 
~dt 



KH 



(14) 



with a natural time-scale tkh = *^sos/— I and an insta- 

' <it Ikh 

bility condition that comes from the classical perturbative 
approach (e.g.. IChandrasekhailll961|) 



PHICM/ 



Phicm/ 



PISM 



PISM 



(15) 



with k wavenumber of the dominant wavelength gas abla- 
tion for Kelvin- Helmholtz instability assumed to be of the 
order of the scale parameter Vs (e.g., iMurrav et aLlll99"3h 
and g = ||g|| from Eqn. [T^l 



3. Synthetic colour magnitude diagrams 

In order to explore the connection between the SFH of 
dwarf galaxies and their physical causes which may be ei- 
ther internal (i.e. star formation itself and feedback) or ex- 
ternal (environment processes: i.e. orbits, tidal field and 
gas dissipative phenomena), we have to deal with the con- 
straints imposed by the observations. The most direct de- 
tection of the SFH is based on the study of CMDs with the 
aid of the synthetic CMDs technique. In brief, given a star 
formation rate (SFR), an initial mass function (IMF), a bi- 
nary fraction and a chemical enrichment law (expressed by 
the helium to heavy element ratio AY/AZ), one can gener- 
ate synthetic CMDs in any photometric system to be com- 
pared with their observational counterparts. To this aim, 
we make use o f the Padova synthetic CMD generator (e.g., 
iNg et all [2002! ). Originally this method was developed for 
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testin g stellar evolution models of stars of dif ferent mass 
(e.g., iBertelli et al.lll98 5l) or stellar clusters (|Chiosi et al.l 
llflSaiBertelh et al.l (l995). We use here the most recent ver- 
sion of this to ol that utiliz es of the library of stellar models 
calculated bv IBertelli et a l. (2009). We refer to the rich lit- 
erature developed by the Padova's astrophysics group for 
further details on YZVAR, the implementation techniques 
and the stellar models in use. 

In the lower right panels of Figs. |31 jSJ IHl we show syn- 
thetic CMDs in the V-I pass-bands without correction for 
the distance modulus. To study the effects of different or- 
bits on the SFH and CMD in turn, we adopt always the 
same model for the dwarf galaxy, evolved along orbits with 
different eccentricities and leave all other parameters un- 
changed. 

By comparing the CMD for the models (E025, E050, 
and E075) it is evident that SFH and associated CMDs 
depend on the orbit eccentricity. The CMD of the E025 



model shows continuity between the episodes of star for- 
mation, which correspond to continuous gas consumption. 
This correlates with the orbital energy of the dwarf galaxy 
in two different ways. On one hand, we have a dependence 
on the location of the dwarf, i.e. on the configuration space. 
The pericentre of the dwarf galaxy is less close to the host 
galaxy than for the other simulations with higher eccentric- 
ity (cf. Fig. [5] and [5]) . According to Eqn. (fTU|) this reduces 
the influence of the external tidal field and coronal gas via 
the potential $ and density p. On the other hand, there 
is a quadratic dependence on the velocity oc , with a 
maximum at the pericentre passage that enhances the ram 
pressure stripping. Therefore, the MW dwarf galaxies with 
smaller eccentricity (say up to e = 0.25) are expected not 
to show clear signatures in the gas consumption of effects 
of dynamical nature in their CMDs. 

Different results are obtained for higher eccentricities, 
cf. the CMD of Fig [51 In this case, we see different subgiant 
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branches as signatures of an irregular consumption of the 
gas. Up to three to four episodes of SF are evident. The 
old metal-poor stars (cyan colour code and small size sym- 
bols) are clearly separated from from the young and metal 
rich ones (red colour code and larger symbols) . Despite the 
overlap of the main sequences (MSs) that should be anal- 
ysed with finer resolution to resolve the sub-components, 
the three to four sub-giant branches (SGBs) are clearly 
indicating three to four episodes. The lowest and faintest 
SGB corresponds to stars of about O.SMq and 9 Gyr among 
which older and lighter stars are also present. This is fol- 
lowed by several episodes, up to ages of about ^ 3.5/4.0 Gyr 
and masses of about m ~ 1.0 AIq. Episodes younger than 
this are difficult to disentangle. It is beyond the aims of this 
study to analyse the CMD in a very detailed fashion. Our 
main purpose here is to show that dynamical effects can in 
principle leave traces of their occurrence in the CMDs of 
the st ellar content of a dwarf galaxy, see [c7^ evic et al.l 
(|201lh and Section [5] below) . 

4. Macrophysics vs. microphysics. An application 
to the Local Group dwarf galaxy Carina 

In the literature there are different approaches to the 
study of the dissipative phenomena in the astrophysics 
context. The motivation of these studies is the attempt 
to explain observational evidence with perturbative ap- 
proa ches to the Navier-Stokes equations (for a review see, 
e.g.. IChandr asckhar 1961). In relation to the astrophysics 
of galaxies, these works have been developed in connec- 
tion with integrated properties (e.g., total mass, central 
velocity dispersion etc.) and s tar formation processes (e.g., 
iLivio et al.lfl980l:lNulsenlll982D . Recently, thanks to the con- 
tinuous development of the numerical techniques, we see 
great efforts in the study of these astrop hysical processes 
from the microsc opic point of view (e.g., iMori fc BurkertI 
l2000l:lQuilisll2004l^ . 

We present in Appendix A a few tests on the compat- 
ibility of the results obtained with our macrospcopic ap- 
proach based on integrated equations and the detailed de- 
scription we can obtain from codes available in the litera- 
ture that solve the microphysics of the instability processes. 
Hereafter, we show the utility of our technique in a real as- 
trophysical context, because the ultimate achievement of 
this technique is indeed to understand the relation between 
the dissipative phenomena and their connection with ob- 
servable quantities, the orbits of dwarf galaxies and their 
CMD. 

The most natural laboratory where to search for real 
observational data is surely the LG of galaxies. The LG 
of galaxies are of paramount importance because of their 
proximity and the wealth of data available so that studies of 
the past SFH are possible fe.g.. iGrebell 119971 : lOrban et al.l 
12008) . In the recent past, the acquisition of large samples 
of spectroscopic data made it possible to measure the in- 
ternal kinematics of some of these dwarf galaxies as well as, 
in a number of cases, their stellar mctallicities and abun- 
dances (e.g.. iKoch et al.ll2008l l2006l 120071) thus initiating 
a new generation of full chemo-dynamical models able to 
take both orbital parameters deduced from proper motions 
and internal kinematics into account. E.g., concerning the 
Magellanic Clouds there is a rich litera ture on t he determi- 
nation of their pro per motions[K allivav alil et al.l (e.g. l2006l) . 
In a recent paper, iPasetto et al.i (|201l1) presented a study 



of the Carina dwarf spheroidal galaxy in which, taking into 
account most of the available observational constraints, we 
recovered in a self-consistent fashion both the past SFH and 
the internal kinematics o f this dwarf ga l axy. 

We use the results of IPasetto et al.l ()2011|) in the light 
of the new technique we have here developed. In this pa- 
per we recover for the first time self consistently the star 
formation history of Carina dwarf galaxy as consequence 
of the pericenter passages which trigger SF. The technique 
that we adopted is based on a chemo-dynamical code (e.g., 
iBerczik et all 120031 ) with new star form ation recipes de- 
scribed in detail in ' Pasetto et al.l ^2010). Our simulations 
reach a mass resolution of 118Mq resolving the core of the 
dwarf galaxy with about 127000 particles (stars) over about 
700000 total particles. In this approach the bursts of star 
formation in the evolution of Carina are the result of micro- 
scopic physical effects: the orbits of particles are integrated 
and some recipes rule the star formation and feedback that 
each particle experiences or exerts on the others. We show 
here that our earlier solution for the star formation history 
of Carina is not unique and it can be easily achieved by our 
new methodology. We p roceed with a few sim pler assump- 
tions than in the work of lPasetto et al.l ()2011f) , because this 
exercise is aimed to be a qualitative example and not a 
quantitative study of the importance of the relative dissi- 
pative effects (that is left for a forthcoming work, Pasetto 
et al 2012 in preparation). We do not spa n the e ntire range 
of orbital parameters as in IPasetto et all (|201lD . moreover 
the starting model adopt ed for the primordial dwarf galaxy 
is inspired by the results I Pasetto et al.l ()2011l ). except that 
for the sake of simplicity we neglect here the temporal evo- 
lution of the model for the MW. The MW potential is de- 
scribed by the following gravitational potential 



MW 
disk 



(i?,z) = - 



GMdisk 



fc4 + v^^+fci 



MW 
bulge 



{R,z) 



bulge 



with G being the gravitational constant and the family of 
parameters {fc^} , i = 2, .., 6, given by {12, 0.8, 6.5, 0.26, 0.7} 
in kpc, and ki = 130.8 in km s'S Mbuigo = 3.4 x W°Mq 
and Mdisk = 1.0 x IQ^^Mq. Over this external field only the 
family of orbits in the neighbourhood of the phase-space 
point X = {22, 89, -34} kpc, v = {-136, 13, -45} in a ref- 
erence frame Sq coUinear with the MW potential are inves- 
tigate and extracted f r om th e family of proper motions in 
Table 1 of lMetz et al.l (2008^ with the minimization action 
technique developed in Paset to et al.l ((2011). Wc adopt an 
initial potential-density pair for Carina from the subfamily- 
(1,4, 7) of the Zhao-models (jZhaolll996h . widely known also 

as the 7- models. We assume M (r) = Mbar (^ip;^) v^iih 

A/bar ~ 4.8 X lO^M© being the baryonic mass of the dwarf 
galaxy, — 0.7 kpc and 7 — 3/2. Finally, we adopt an 
initial SFR Vo = O.O29M02/r~^ at ^ib = lOGyr (iib stands 
for look-back time). The corresponding initial value of the 
pressure required to balance the cloud destruction rate and 
the star formation rate is P = O.OIP©. The results for the 
recovered SFH are plotted in Fig. [71 Finally, stellar mod- 
els in system Eqn. ^ and for the H-R diagrams in the 



10 



S. Pasetto et al.: Dissipative phenomena in LG dwarf galaxies 



QC 3 

Ll. 

w 

CO ^ 




-6 -4 

t„ 



-2 



Fig. 7. Recovered episodic star formation history of the 
Carina dwarf galaxy superposed_ on the normahzed his- 
togram bars from th e data of iRizzi et alJ (|2003h (see also 
iPasetto et all (|2Qllf )V 



low mas s regime (from 0.15 to 2.5 Mq) come from the 
work of BcrtcUi et alJ ([2003) , and for higher mass from 
iBertelli et ah (20091). The technique is based on the the e ar- 
lier works oflBressan et alJ (|1993D . lFagotto et all ()1994allb[ ). 
iGirardi et alT ( T996^ but is updated in the opacities, rates 
of energy loss by plasma neutrino s and nuclear network in- 
tegration f see IBertelli et aLll2009l and references therein). 

By straight c o mpar ison of our plot in Fig. [7] with Fig. 11 
in iPasetto et al.l (|201lf ) , the similarity is remarkable, thus 
confirming the original idea contained in that paper that 
the SFH of the Carina dwarf galaxy can be explained as 
results of pericentre passages. But the techniques are ab- 
solutely different: In our approach we do not integrate the 
orbits of the single particles. There are not microscopic par- 
ticles but macroscopic analytical properties (total mass, 
density profiles, etc.). The star formation is the result of 
the pressure on the molecular cloud model, and the pres- 
sure depends on the orbit thought the MW potential, the 
tidal tensor as well as ram pressure, Kelvin-Helmholtz and 
Ray leigh- Taylor dissipative phenomena. The major differ- 
ence is that the effect of the tidal field of the MW is fully ac- 
counted for by solving the EoM for the single dwarf galaxy 
particles in a N-body simulation, while here it is considered 
only by means of a first order perturbation, ultimately by 
the tidal tensor T. More importance is given here to exter- 
na l instability processe s that were not taken into account 
in IPasetto et al.l (|201lh . We have two different techniques 
that match the same observational constraint, i.e. the star 
formation history. We point out another im portant differ- 
ence w ith respect to the results obtained by IPasetto et al.l 
()2011|) : while a high resolution chemo-dynamical simula- 
tion requires about 4 months on a special hardware ma- 
chine, the result presented here needed only about 8 hours 
on a parallel cluster of processors. Moreover, the present 
technique can impose more straightforward constraints on 
the gas fraction adopted for modelling Carina that are im- 
possible to reach because of particle resolution problems 
in a fully hydrodynamical NB-TSPH simulation. In this 
respect the new tec hnique is superior to that followed by 
IPasetto et al.l (|2011[ ). The high speed of the present algo- 
rithm is partially due to the fully analytical formulation 



and partially to the high scalability of the code. We refer 
the interested reader to Appendix A for details on the par- 
allelization performance of the code. Alternatively, the fully 
hydrodynamic code is best suited to studies of the internal 
kinematics that we have not taken into consideration here. 

5. Discussion, open questions, future developments 
and conclusion 

We have developed a method and companion code able to 
handle in a fast way the major phenomena affecting the 
SFH of extended dwarf galaxies moving through an ex- 
ternal environment and compared the resul ts with those 
obtained from a fully hydrodynamical code (jPasetto et al.l 
I2OIII ). This formalism stands on a novel interpretation of 
the classic al results for a flow past an extended spher ical ob- 
ject ie.s... lBatcheIo3l2000l: iLandau fc Lifshitdll959D where 
we considered the problem in a non-inertial reference frame 
immersed in an external field and we exploited the tidal ap- 
pr oximation for the external gravitational interaction (e.g., 
iPasetto fc Chiosi 2009. ) . 

The advantage of this method, if compared to fully hy- 
drodynamic codes, is its high speed (without imposing the 
restriction of a point-mass approximation) which is possi- 
ble thanks to a new direction-dependent formulation for the 
pressure effects described in a non inertial reference system 
and that is a valuab le a lternative to t he classical works of 
iGunn fc GottI ([l97l or lNulsenI (|T98l . 

The technique is very suitable to impose constraints 
on the orbits based on the SFH of dwarf galaxies. We 
have shown a possible application to the case of Carina 
to confirm the validity of the approach adopted. Another 
natural applications is the study of the triplet Large 
Magellanic Cloud, Small Magellanic Cloud and MW by tak- 
ing into account their SFH. A further possible application 
is to constrain the phase-space parameters of the external 
group of dwarf galaxies. The apparent lack of intermedi- 
ate age populations in some early-dwarfs around Cen A 
(see Crnoicvic ct al. 2011,) or in other objects for which, 
due to their distance, only the most luminous part of the 
CMD can be resolved into stars could be investigated using 
the structural and dynamical constraints we have proposed 
here. In studies of such more distant dwarfs, spectroscopic 
measurements of internal kinematics are still difficult to 
obtain and often only photometric data are available. As 
a consequence, the investigation of these important issues 
with complete chemo-dynamical codes is time-consuming 
and time consuming, whereas a fast and flexible technique 
and companion codes have better chances. 

Finally, while the competition of the different dissipa- 
tive phenomena to produce an observable quantity is surely 
a target of this project, it is not the primary goal. A dwarf 
galaxy is affected in its orbital evolution by several phe- 
nomena that all act jointly to produce an observational 
result, making it difficult to constrain the contribution of 
one individual mechanism. Our technique is, by itself, not 
able to disentangle the relevance of all those different phe- 
nomena because, as just seen in the exercise for Carina 
dwarf galaxy, the same result can be achieved with differ- 
ent methods. Nevertheless, this approach is a step forward 
in the inclusion of the stellar population algorithms to add 
a further observational constraint to the orbit investigation. 

Of course, the simplicity of the adopted approach comes 
also with a few technical limitations and drawbacks that 
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should be amended. For instance, in the Ra yleigh- Taylor in- 
stability the magnetic field can play a role (IChandrasekhail 
[l961). Moreover the magnetic fields might also affect the 
star formation efficiency by changing the Jeans instability 
criteria because the critical mass for collapse depends on the 
Alfven and the sound speeds. Even if magnetic fields play 
an important role in the collapse and fragmentation of cold 
clouds, the generalization of Eqn. pU)) to include magnetic 
fields is not a trivial task. Moreover, if we consider that the 
interstellar medium in dwarf galaxies is heterogeneous and 
has varying density, we expect hydrostatic turbulent pres- 
sure to act in differently w ay in different regions of a dwarf 
galaxy ([Bureau fc Carign"an ,2002' ). Finally, tightly related 
to the studies of star formation and environmental effects is 
the problem of gas removal from dwarf galaxies. If ram pres- 
sure and/or gas consumption cannot explain all the present- 
day observations, a plausible alternative could be the gas 
ionization by UV radiation (of stellar or Galactic origin) 
and by SN explosions. These effects played an important 
role during the epoch of reionization by increasing the the 
Jeans mass. Unfortunately the theoretical predictions for 
the low mass regime of the dwarf galaxies is strongly bi- 
ased by the numerical resolution. Finally, in a cosmological 
context the reionization at low redshift can hav e played an 
impo rtant role in the star formation processes (jlliev et al.l 
|2011| ). which however we have neglected here to first ap- 
proximation of the problem. 

Considering our finding, a natural follow-on project 
would be a study of an extended family of orbits and model 
parameters to investigated the competitive role of the dif- 
ferent dissipative phenomena acting around the MW, M31 
or by considering the LG tidal interactions: which processes 
affect the star formation history, and which do not? When 
is a deep potential well of the dwarf galaxy sufficiently mas- 
sive to suppress the effects of the various dissipative phe- 
nomena? When do tides promote or inhibit star formation 
through compression/stripping? Similarly for ram pressure, 
Ray leigh- Taylor, Kelvin-Helmholtz effects, when and where 
are they dominant in the history of MW dwarf galaxies? 
These and many other questions can be investigated with 
this technique, for the first time in relation with syntheti- 
cally generated colour magnitude diagrams, in order to un- 
derstand how to interpret the observed CMDs (Pasetto et 
al. 2012, in preparation). 
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Appendix A: Testing the code 

We group in this section a few technical tests performed to 
control the code stability, its convergence and to compare 
with other literature results. 



A.l. Convergence tests and code scalability 

We present here the test on the convergence of the re- 
sults of the integration of the system of equations proposed 
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Fig. A.l. Plot of theoretical load balance between number 
of processors Np and number of equations in the system 



in Eqns. ^ and a few comments on a remarkable result 
we have obtained with the parallelization of the code. As 
explained in the main text, our problem is to single out 
the volume occupied by the orbits (differentiable functions 
a; : K ^ M"^ of the real line (i.e., the time) on the con- 
figuration space Q) that minimizes a function-difference 
between an observed and a synthetic CMD (convolved with 
distance modulus and errors). This function is clearly sur- 
jective (i.e. 'onto') so that, once the volume of the param- 
eter space is huge and multidimensional, the possibility of 
rapidly converging to a best-fit solution becomes an impor- 
tant issue. The structure of the code we have developed is 
extremely promising thanks to the scalability of the sys- 
tem of Eqn. For its integration the classical Runge- 
Kutt a method w i th ad aptive step-size has been adopted 
(e.g., iPress et all Il993l) . but the mass classes have been 
spread over a large number of processors. Clearly to op- 
timize/balance the load we need to split the equations over 
the rank size by minimizing the latency over the available 

processors: — Int 



Np with Np the number of pro- 
cessors and Nc the number of classes. The plot is as in the 

Fig.im 

Once the performance of the code was understood, we 
performed a few convergence tests on the astrophysical re- 
sults. Depending on the available computational resources, 
say the Np, we can choose the best number of processors 
on which to perform our integration by looking at Fig. lA.ll 
We had the possibility to test the lower minimum peak on 
the right of Fig. [XT] (for Np = 6500). For this high reso- 
lution test, the system of equations of Eqn. (U) has been 
integrated with a ~ as defined in Section 12.11 The 

results are presented in Fig. IA.21 where four lines are plot- 
ted, for the integration of the system of equation Q with 7, 
65, 650, 6500 molecular cloud classes on an E050 orbit. The 
test shows how the approximation adopted in our study, for 
a = 1/100, i.e., 650 molecular classes (the blue line), gives 
a robust consistent star formation rate with the higher res- 
olution simulation (the red line): red and blue line overlap 
almost perfectly. 

On 1st May 2011 the simulation ran on 6500 processors 
at SGI Altix 4700 at the National Supercomputer HLRB- 
II, Munich ( Germany) on a dedicated queue recording a 
peak performance of 23Tflops/second. 
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A.2. Tests on the integro-difFerential system of equations 
Eqn. (HI) 

Eqn. is an integro-dilFerential system of equations 
used to deduce the star formation rate of a galaxy from 
the total pressure P. While this system is conceptually 
stable because it is based on the funda mental princi- 
ples o f stel lar structure and evolution (e.g., iBertelli et al.l 
[2001 [2OOI . its integration is not trivial and the numer- 
ical approaches adopted needs to be tested. Its integra- 
tion on the independent variable t is the most computa- 
tional time consuming part of the code and has been per- 
formed with distributed-memory parallelization standard 
techniques (message-passing interface, MPI). A serial in- 
tegration of a si milar system has been already performed 
in literature by iFuiital (|l998l) . Therefore, we can p repare 
a case -test to reproduce the same result posted in iFuiital 
()1998[) by omitting the pressure determination with our 
Eqn. ([TOl) and instead by "injecting" in our code t he syn- 
thetic pressure profile of Fig. la (dashed line) in IFuiital 
(fT99l . Thus, we eliminate the code section devoted to our 
original pressure determination and only the parallel inte- 
grator of the system of equation Eqn. (jH) is tested against 
the ser i al ind ependent determination originally presented in 
IFuiital (|1998ll . The pressure is defined piece- wise as follow: 



Pit) 



Pq t<0 

p^gf(ln2+ln5)i Q < t < 2 

IOOOPq t > 2. 



(A.l) 



We see in Fig |A.3l that our results that are almost identical 
to the dashed line of Figure lb in IFuiital (|l998) . Here the 
integration is performed with Runge-Kutta met hod with 
adaptive stepsize control (e.g.. iPress et al.lll993l ) and dis- 
tributed over 650 processors. 



A. 3. Pressure from Ean. (fTOj) and and speed limits 

We assume here that the dwarf galaxy is falling in an gen- 
eral intra cluster medium (ICM) or that a dwarf galaxy 
is passing through the MW gaseous disk. In these situa- 
tions the motion of the dwarf can easily be in supersonic 



regime, i.e. Mach number Mpre > 1, or equivalently that 
the flow is impacting on our dwarf galaxy at supersonic 
velocity Vg > Vsjcm with WsJCM sound speed of the flow. 
Moreover, we keep general this exercise assuming a poly- 
trophic equation of state equation with adiabatic index 7 
for the gas flowing. The presence of a shock increase density 
(and pressure) by compression while the velocity field from 
supersonic become subsonic (we will limit our arguments 
to normal shocks) whence Ean. (|TU|) holds. Nevertheless, it's 
simple to prove that the pressure of the gas before the bow 
shock, Pprc, and the stagnation point Ps can be very differ- 
ent and the Eqn.([TU]) have to be 'clothed' with supersonic 
formalism. We call the pressure after the surface of discon- 
tinuity (the thin shock), Ppost- In this case the relations 
betw een pressures pre and po st shock are already known 
(e.g., iLandau fc Lifshitdll959 l ) 



P 



post 



J pre 



1 



27 



7 + 



1 



7 + l + (7-l)(M2,.,-l) 



(A.2) 



7+l + 27(M2,,-l) 
and we obtain the pressure between Ppost cind Pg as 

p. 



p 



post 



(A.3) 



With M2^,t and 
tion between 



Ppro 

■ as: 

P.. 



from Ean lA.2l we can write the rela- 



Ps P 



post 



P 



pre 



P 



7-1-1 



post Pprc 

27 
7 + 1 



Mi 



_^ (A.4) 

7-1 



whose plot is as in Fig |A.4l As evident, the pre-shock pres- 
sure and the pressure at the stagnation point are not of 
the same order even if the post shock velocity field is sub- 
sonic: the flow is compressed and its pressure (and density) 
increases passing through the shock and after the shock 
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Fig. A. 4. Tend between of the ratio between pressure at 
the stagnation point and pressure pre shock as a function 
of the Mack number of the gas where the dwarf galaxy is 
moving. 



before reaching the stagnation point. As seen in Fig. IA.4I 
these pressure remain of the same order only for subsonic 
or weakly supersonic that represent the limit of validity 
of Ean. (fTU|) . before supersonic correction has to be imple- 
mented. 

Appendix B: Star formation efficiency and 
disruption time relation 

The connection between star formation processes and 
molecular clouds is a ver y fertile research topic encom- 



passing observations (e.g.. BTgiel et al.1 ^2008'). theory (e.g., 
iKrumholz fc McKee 20051) and experimental/numerical 
works (e .g.. 1X168861112011). In our app roach we revisited a 
work of lElmegreen fc Efremovl ()1997[ ) chosen for its sim- 
plicity. Other literature results can be similarly imple- 
mented. Let Mi be the cloud's initial mass as defined in 
Section [2.1l for which we assumed a constant sta r formation 
rate c in the integration interva l dt (see also Elraegreenl 
[1981 IKrumholz & McKed [MqI IKrumholz fc TanI I2007D . 
Following Surdin (198^, the erosion rate is proportional to 
the luminosity Ls of embedded stars with total mass Mgtar.i 
divided by the specific mass binding energy Mio'^ with a 
dispersion velocity. The equation for the rate of change of 
gas mass in the cloud considered is: 



the luminosity 



dMi.__ _ A ^ 

dt ~ V2' 



U {t) = / <;^A^ {t ~ t') dt' 

Jo 



(B.l) 



(B.2) 



where (t) is the luminosity-to-mass ratio of a population 
of stars generated by the cloud class at the instant t that 
we obtained as explained in Section |21 Let us drop for the 
moment the subscript i to simplify the notation. From the 
previous Eqns. (jB.ll) and (IB.2[) we get: 



dt' / A{t'- t") dt" 



which can be integrated to give 
M (t) = M (0) - ? ( t + 4 



where the double integral is more easily numerically com- 
puted as: 



to give: 

M{t) 



f dt' f A{t'-t")dt" = 
Jo Jo 



A{t~t')dt'- / t'A{t')dt' 



A - t') dt' 



(B.5) 



M(0)-df + 



-^*t'A {t')dt'^^ 



(B.6) 



for every i. We can find also a quicker approach based on an 
interpolation function in the work of Girardi et al.l (Il995l 
their Fig. 13) as recently in Elmcgrce n et al.l (|2006[ ) that 
we can also adopt when necessary to speed up the code. 

By using a power-law like A (i) = Aq ^-^^ , A e [0, 1[, the 

integrals in Eqn. (|B.6p can be carried out analytically. We 
observe that the destruction time r of the molecular clouds 
can be defined in an implicit way as the instant i\M (t) = 0: 



At^Aoc; 



(A-1)(A~2), 



— -I- M (0) = ^ 



1 



ar 



-2-A 



0, 



(A-l)(A-2) 
where we exploited the following definitions: 

t 



and 



to^ 



M(0) 
^Ao M (0) 



(B.7) 



(B.9) 



with to ~ lO'^ yr as a fixed parameter from I Girardi et all 
(|1995D and A ~ 0.6 for t > to and A = for t < to where the 
solution is analytical (a quadratic equation). Apart from 
the difference in the mathematical formalism here laid out, 
the co ntents follow then exactly as in Elmegrcen & Efrcmc^ 
(|1997D to which we refer the reader for extended discussion. 
We point out here that we can define the efficiency as 

The two functions defined in (IB. 81) are of interest due 
to their dependence on pressure impacting on the or- 
bital mass of the dwarf galaxies in o ur LG. We as- 
su me the following functiona l dependence (jElmegreenll 1 9891 : 
El megreen fc Efremovlll997l) : 



= ao I 



a = a (M, P) 



M 



ViOHfo 
M 



-1/4 



-1/4 



p 
p 



3/8 



-5/8 



(B.ll) 
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with dimensionless constants ao = 0.1, /3o — 180, Pq = 
3 X 10'^A;Bcm~^K, and ks as Boltzmann constant. 
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